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ABSTRACT 

Solar flare X-ray emission results from rapidly increasing temperatures and emission measures in 
flaring active region loops. To date, observations from the X-Ray Sensor (XRS) onboard the Geo- 
stationary Operational Environmental Satellite ( GOES) have been used to derive these properties, 
but have been limited by a number of factors, including the lack of a consistent background sub- 
traction method capable of being automatically applied to large numbers of flares. In this paper, 
we describe an automated temperature an d emission measu re-based background subtraction method 
(TEBBS), which builds on the methods of lBornmannl ([19901 ). Our algorithm ensures that the derived 
temperature is always greater than the instrumental limit and the pre-flare background temperature, 
and that the temperature and emission measure are increasing during the flare rise phase. Addi- 
tionally, TEBBS u tilizes the improved estimates of GOES temperatures and emission measures from 
IWhite et al.l (12005). TEBBS was successfully applied to over 50,000 solar flares occurring over nearly 
three solar cycles (1980-2007), and used to create an extensive catalog of the solar flare thermal prop- 
erties. We confirm that the peak emission measure and total radiative losses scale with background 
subtracted GOES X-ray flux as power-laws, while the peak temperature scales logarithmically. As 
expected, the peak emission measure shows an increasing trend with peak temperature, although the 
total radiative losses do not. While these results are comparable to previous studies, we find that 
flares of a given GOES class have lower peak temperatures and higher peak emission measures than 
previously reported. The resulting TEBBS database of thermal flare plasma properties is publicly 
available on Solar Monitor (www.solarmonitor.org/TEBBS/) and will be available on Heliophysics 
Integrated Observatory (www.helio-vo.eu). 



1. INTRODUCTION 

Solar flares are among the most powerful events in the 
solar system, releasing up to 10 33 ergs in a few hours or 
even minutes. They are believed to be powered by mag- 
netic reconnection, a process whereby energy stored in 
coronal magnetic fields is su ddenly released. According 
to the CSHKP flare model (iCarmichaeJ H96l iSturrockl 
[19661 [ffirayama][l97l IKodd fc P neumanll 1 9 761 ) . electrons 
accelerated by magnetic reconnection spiral down the 
magnetic loops and strike the chromosphere causing the 
emission of hard X-rays (HXR). As a consequence, the 
chromospheric material is also heated and expands back 
up into the loops which causes the obser ved increase in 
temp erature and emission measure (e.g., iFletcher et al.l 

To date, the study of solar flares has been predom- 
inantly focused on single events or small samples of 
events. While such studies have furthered our under- 
standing of the physics of these particular flares, they 
are fundamentally limited since they cannot, with any 
certainty, explain the global behavior of solar flares. In 
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contrast, only the study of large-scale samples can give an 
insight as to whether findings of given studies are partic- 
ular to individual events or characteristic of many. This 
can allow constraints to be placed on global flare proper- 
ties and give a greater understanding to the fundamental 
processes which drive these explosive phenomena. 

That said, large-scale studies of solar flare properties 
have been few in number o ver the past decades . Such 
a study was performed by iGarcia fc Mcintosh] (1992) 
who used the Geostationary Operational Environmental 
Satellite (GOES) to examine 710 M- and X-class flares. 
They noted a sharp linear lower bound in the relation- 
ship between emission measure and GOES class. How- 
ever, this paper is mainly focused on categorizing types 
of very high temperature flares and examined whether 
these flares approached or exceeded this emission mea- 
sure lower bound. 

A definitive example of a large-scale study of 
the thermal prope rt ies of solar flares was conducted 
by iFeldman et al.l (|1996bl ). who combined results 
from three previous studie s ([Phillips fc Feldmanl li~995l: 
IFeldman et al.l 119951 Il996al ) to investigate how temper- 
ature and emission measure vary with respect to GOES 
class for 868 flares, from A2 to X2. Their work used tem- 
peratures derived using the Bragg Crystal Spectrometer 
(BCS) onboard Yohkoh. These temperature values were 
convolved with the corresponding GOES data to derive 
values of emission measure. They found a logarithmic re- 
lationship between GOES class and temperature, and a 
power-law relationship between GOES class and emission 
measure, with larger flares exhibiting higher tempera- 
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tures and emission measures. However, temperature and 
emission measure were derived at the time of the peak 
1-8 A flux and so are likely to be less than their true max- 
ima. Furthermore, BCS temperatures hav e been found to 
be hig her than those measured by GOES (Fcldman et al. 
1996b), and using these values to calculate GOES emis- 
sion measure will give lower values than if GOES was 
used consistently 

More recently, iBattaglia et al.l (|2005| ) studied the cor- 
relation between temperature and GOES class for a sam- 
ple of 85 flares, ranging from Bl to M6 class. Al- 
thou gh the values r eporte d gave a flatter dependence 
than lFeldman et all (|1996bl ) , the large scatter in the data 
led to a very large uncertaint y making the two rela - 
tions comparab l e. In c ontrast to lFeldman et al.1 |T996b), 
IBattaglia et al.1 (|2005l ) accounted for solar background 
and extracted the flare temperature at the time of the 
HXR burst as measured by the Ram aty High-Energy 
Solar Spectroscopic Imager (RHESSI; ILin et all I2002f) 
rather than at the time of the soft X-ray (SXR) peak. 
However, any discrepancies expected to be caused by 
these differences were not di scernible in y iew of the large 
uncertainties. In addition, iCaspil (|2010fl used RHESSI 
to examine the temperature of 37 high temperature 
flares with GOES class. The relationshi p was quali- 
tatively similar to tho se of iFeldman et al.1 (|1996bl ) and 
IBattaglia et al.1 (|2005h however as the relationship was 
not fit no quantitative comparison can be made. 

Larger st a tistica l sa mples were studie d by 
iChriste et all (|2008ft and lHannah et all ([20081 ) . who 
investigated the frequency distributions and energetics 
of 25,705 microflares (GOES class A-C) observed by 
RHESSI hom 2002 to 2007. From those events for which 
an adequate background subtraction could be performed 
(6,740) a median temperature of ^13 MK and emission 
measu re of 3xl0 46 cm -3 were found. lHannah et al.l 
(2008), in particular, looked at the temperature derived 
from RHESSI observations as a function of (background 
subtracted) GOES class, and foun d si milar trends to 
the wo rks of IFeldman et al.1 (|1996bT ) andfBattaglia et al.1 
(2005). However, their analysis only included events of 
low C-class and below. 

While these studies have provided some insight into 
the global properties of solar flares, they each have their 
limitations. In particular they lack a commonly used 
method of isolating the flare signal from the solar and 
instrumental background contributions. Previous back- 
ground subtraction methods have often been performed 
manually. Others, such as setting the background to 
the flare's initial flux values, or fitting polynomials be- 
tween the flux values at the start and end of the flare, 
often exaggerate noise and do not preserve characteristic 
temperature and emission measure evolution. Therefore, 
the accurate separation of flare signal and background 
limits the number of events tha t can be analyzed. For 
example, IBattaglia et al.l ([2005), in accounting for so- 
lar background, were only able to compile a sample of 
85 events. Although a larger dataset would not have 
reduced the range of scatter, it would have better re- 
vealed the variations in the density of points within the 
distribution. This would have allowed a fit to be more 
tightly constrain ed and thereby reduced the uncertain- 
ties. Conversely. IFeldman et all (|1996b|) . with a sample 



of hundreds of flares, did not attempt to account for the 
solar background at all, which can bias smaller events as 
the background makes up a greater contribution to the 
overall flux. 

Few attempts have been made to develop automated 
background subtraction techniques for GOES observa- 
tions which can b e applied to large numbers of flares. 
IBornmannl (|1990h developed a method to determine 
whether a given background subtraction preserves char- 
acteristic temperature and emission measure evolution 
without checking manually, i.e., that temperature and 
emission measure both increase during the rise phase 
of flares. This beh aviour has been seen in numer- 
ous o bservations (e.g.. lFludra et al.l [l995; Battagli a et all 
20091) and numerical models (e.g.. IFisher et a~ 119851 : 
Aschwanden fc Tsiklauril 120091 ) . Th is method used the 
polynomials of iThomas et all ([1985) which relate tem- 
perature and emission measure to the ratio of the short 
and long GOES channels, R = Fs/Fl- However, 
IWhite et all (H005) have since improved on this by as- 
suming more modern spectral models (CHIANTI 4.2 
ILandi et al.11199 9. 2002) and taking into account the dif- 
ferences between cor onal and photosph eric abundances, 
requiring the tests of IBornmannl (|1990f ) to be updated. 

In this paper, we study the thermal properties of solar 
flares using GOES observations over nearly three solar 
cycles. The flare signal within these observations has 
been isolated from the various solar, non-solar, and in- 
strumental background contributions using a modified 
background subtraction method. This in turn has al- 
lowed more accurate automatic calculation of flare prop- 
erties. In Section [5] of this paper, we discuss the GOES 
XRS, the GOES event list and how to derive plasma 
properties from GOES observations. In Section [3] we 
describe previous background subtraction methods for 
GOES observatio ns and out l ine ho w we have improved 
upon the work of IBornmannl (|1990l ). In Section 0] we use 
this method to improve upon previous statistical stud- 
ies by deriving flare properties such as peak temperature 
and emission measure for flares in the GOES event list 
and examining the relationships between them. In Sec- 
tions [5] and [6] we discuss the results and provide some 
conclusions. 

2. OBSERVATIONS & DATA ANALYSIS 

2.1. Geostationary Operational Environmental 
Satellite /X-Ray Sensor 

The observations used in this study have been made by 
satellites of the GOES (Geostationary Operational Envi- 
ronmental Satellite) series, which has been operated by 
the National Oceanographic and Atmospheric Adminis- 
tration (NOAA) since 1976. Each spacecraft carries an 
X-Ray Sensor (XRS) onboard which measures the spa- 
tially integrated solar X-ray flux in two wavelength bands 
(long; 1-8 A, and short; 0.5-4 A) every three seconds. 
The sensitivities of the various XRS instruments have 
remained comparable over the years, although the de- 
sign of the GOES-8 XRS and subsequent detectors was 
altered due to the change from a spin-stabilized to 3-axis 
stabilized platfor m. For an in-dep t h disc ussion of the 
GOES-8 XRS see IHanser fc Sellers! (fl996h . The GOES 
series has provided a near-uninterrupted catalog of so- 
lar activity for over three complete solar cycles, and the 
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Fig. 1. — X-ray lightcurves of an M1.0 solar flare observed by 
GOES, a) X-ray flux in each of the two GOES channels (0.5-4 A; 
dotted curve and 1-8 A; solid curve), b) The derived temperature 
curve, c) The derived emission measure curve. The vertical dotted 
and dashed lines denote the defined start and end times of the 
event, respectively. The vertical red, black and green lines mark 
the times of the peak temperature, peak 1-8 A flux, and peak 
emission measure, respectively. 



GOES flare classification scheme is now universally ac- 
cepted. 

2.2. The GOES Event List 

The sample of flares used in this study has been ex- 
tracted from the GOES event list; a list of solar X-ray 
events which has been compiled by NOAA throughout 
the lifetime of the GOES series. In order for a solar flare 
to be included in the GOES event list, it must satisfy two 
criterigQ: firstly, there must be a continuous increase in 
the one-minute averaged X-ray flux in the long channel 
for the first four minutes of the event; secondly the flux 
in the fourth minute must be at least 1.4 times the initial 
flux. The start time of the event is defined as the first 
of these four minutes. The peak time is when the long 
channel flux reaches a maximum and the end of an event 
is defined as the time when the long channel flux reaches 
a level halfway between the peak value and that at the 
start of the flare. 

The flare start and end times determined by these def- 
initions do not always agree with those identified manu- 
ally. An example of this can be seen in Figure [TJa which 
shows the X-ray fluxes in the two GOES channels for 
an M1.0 solar flare that occurred on 2007 June 2. The 
event list start and end times are marked by the vertical 
dotted and dashed lines, respectively. The start time of 
the GOES event is a couple of minutes before the on- 
set of the flare. Nonetheless, this start time satisfies the 
event list criteria and highlights a drawback in the event 
list definitions. Another drawback is associated with the 

7 http: //www. swpc .noaa.gov/ftpdir/indices/events/README 



event list end time. It can clearly be seen that the decay 
of the flare in Figure QJi continues for over half an hour 
after the event list end time. This means that proper- 
ties depending on the decay time or duration of the flare, 
such as total radiative losses, will be systematically un- 
derestimated. 

However, these definitions also help reduce the number 
of 'double-flares' in the event list, i.e., two flares being 
incorrectly labeled as one. This can happen when one 
flare occurs on the decay of another thereby preventing 
the full-disk integrated flux reaching half the peak value 
of the first flare. Having searched the event list between 
1991 and 2007 we found 1,865 out of 34,361 events (5.4%) 
contained points between their peak and end times which 
satisfied the event list start criteria. Of these, the second 
flare was recorded in the event list in 236 cases. It should 
also be stated that the event list start criteria do not 
locate small events (e.g., B-class) at times of high back- 
ground flux or during large flares (e.g., M-class). This 
is because a small flare will not cause the full-disk inte- 
grated X-ray flux to increase to 1.4 times the initial value 
when that initial value is more than an order of magni- 
tude greater than the flare itself. Therefore, although 
one would expect to always find more small events, the 
event list actually contains fewer around solar maximum 
when large events are more frequent and the background 
is often at the CI level or higher. 

The GOES event list for the period 1980 to 2007 was 
used in this study. Data from the 1970s were not included 
due to their poor quality and because many GOES events 
from this period were erroneously tagged. This meant 
that a total of 60,424 events, from B-class to X-class, 
were considered. Events for which data were unavail- 
able, erroneously included (i.e., did not satisfy the event 
list definitions), or displayed data drop-outs were then 
removed. It was found that the size distribution of the 
discarded events was very similar to that of the entire 
data set. This implies that the remaining dataset was 
not biased by the exclusion of these events. After these 
events were removed, 52,573 remained. 

2.3. Deriving Flare Plasma Parameters 

Although GOES only measures X-ray flux in two pass- 
bands, techniques have been developed to derive flare 
plasma prope rties from the ratio of the short a nd long 
channels (e.g. iThomas et alJll985t lGarcialll994[ ). These 
properties include temperature, T, emission measure, 
EM, and the total radiative loss rate from the X-ray 
emitting plasma, dL ra d/dt. In this study, temperature 
and em ission measure were computed from the calcula- 
tions of I White et al.l (|2005| ). who used updated detector 
responses and plasma source functions to create tables 
of the dependence of temperature and emission measure 
on the fluxes in GOE S channels. This meth od is an up- 
dated version of that of I Thomas et al.l (|1985T ) who derived 
temperature and emission measure relations from poly- 
nomial fits to the GOES-1 XRS response function. The 
tabulate d values of temper ature and emission measure 
given in iWhite et al.l ()2005l ) can be approximated using: 
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where R(= Fs/Fl) is the ratio of the short to long chan- 
nel. For values of the coeffi cients A n and B n for each 
GOES satellite, see Table 2 of lWhite et all (120051) . These 
tables cover the range from 1-100 MK. However, due to 
instrumental sensitivities of the XRS instruments used 
in this study they are only valid above 4 MK. 

Figures QJ> and [T]c show the temperature and emission 
measure evolution of the 2007 June 2 flare. Their be- 
havior is characteristic of the evolutio n of a typical flare . 
According to the standard flare model (|Carmichael | 1964 |: 
ISturrocklll966t IHiravamalll974t IKopp fc Pneumanlll976l 
CSHKP), nonthermal electrons rapidly heat the chro- 
mosphere to high temperatures (Figure QJ>) , causing the 
plasma to expand into overlying flare loops in a process 
called chromospheric evaporation. These high loop den- 
sities result in an increased emission measure following 
the temperature peak (Figure [T]c). Once heating has 
ceased, the plasma cools by thermal conduction and then 
via radiative processes. This is accompanied by a pro- 
gressive decrease in both flare temperature and emission 
measure. 

Theoretically, the radiative loss rate, dL rQ( z/cft, can be 
calculated using 



dL r 



dt 



EM x A(T) erg s" 



(3) 



where A(T) is the radiative loss function. Here, this was 
estimated using tables of radiative loss rate as a function 
of emission measure for various te mperatures, which were 
generated using CHIANTI (v6 llDere et al.ll2009l ) and 
the methods of ICox fc Tuckerl ([19690 . Tni s technique 
may lead to under- or over-estimates of the true radiative 
loss rate because it assumes an isothermal plasma which 
may not well approximate the flare's differential emission 
measure (DEM) distribution. 

In the above ca lculations, coronal abundances 
dFeldman et al. 119921) and the ioniz ation equilibria o f 
Mazzotta et al. (1199 81 were assumed. iDere et al.1 (2009) 



justified the use of these equilibria by compa ring them to 
others obtained from the ionization rates of iDerd ([2007D 
and a revised set of recombination rates. The results were 
found to be similar. A const ant density o f 10 10 cm~ 3 
was also assumed justified by iWhite et al.l (|2005l ) who 
used CHIANTI to compute the spectrum of an isother- 
mal plasma at 10 MK with densities of 10 9 , 10 10 , and 
10 11 cm~ 3 , and found no significant differences between 
them. 

Finally, having calculated the radiative loss rate, the 
total radiative losses of the X-ray emitting plasma in the 
flare can be calculated by integrating between the flare 
start and end times: 



Erad — 



dL ra d(t) 

dt 



dt ergs 



(4) 



3. BACKGROUND SUBTRACTION METHOD 



As the GOES lightcurves do not include any spatial 
information, they contain contributions not only from 
the flare but also from all non-flaring plasma across the 
solar disk. In addition, the lightcurves include non-solar 
contributions such as instrumental affects which vary be- 
tween the individual X-Ray Sensors. These various back- 
ground contributions can cause significant artifacts when 
deriving flare properties. Therefore, it is imperative to 
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Fig. 2. — Schematic of a flare X-ray lightcurve showing how the 
total flux detected by, e.g., the GO ES XRS, is divide d into con- 
stituent components. (Adapted from Bornmann 1990). The total 
flux (solid line) is the sum of the flux from the flare plus the so- 
lar background (divided by the dashed line). The pre-flare flux, 
however, is the sum of the background component and the quies- 
cent component of the flaring plasma (e.g., the associated active 
region). 

isolate the flare signal from these contributions, partic- 
ularly for weaker events. In this study a generic tem- 
perature and emission measure-based background sub- 
traction method (TEBBS) ha s been developed w hich im- 
proves upon the methods of iBornmannl (|1990t ). These 
methods aim to discard possible flare signals which do 
not preserve the increasing nature of temperature and 
emission measure during the flare's rise phase (high- 
lighted in Section [lj . In this section, the limitations of 
previous background subtraction methods are discussed, 
before the TEBBS method is described in detail. 

3.1. Previous Background Subtraction Methods 

The schematic in Figure [2] shows how a hypotheti- 
cal GOES lightcurve is divided into its flaring and non- 
flaring (i.e., background) components. The two limit- 
ing cases in calculating the the boundary between back- 
ground and flare fluxes are either to assume that the to- 
tal flux is dominated by the flare, thereby not performing 
any background subtraction, or to assume that the back- 
ground is equal to the flux near the beginning of the event 
('pre-flare' flux). The first assumption may be valid for 
events which are orders of magnitude above the back- 
ground level, but is clearly incorrect for weaker events. 
The second assumption may be incorrect as there may 
be significant flare emission before the flare detection al- 
gorithm reports the sta rt time. An example o f the first 
method can be found in lFeldman et al.l (|1996b| ) in which 
no background was subtracted. An example similar to 
the second method can be used in the GOES workbench^ 
which allows the background to be calculated as a line 
(polynomial or exponential) between the flux values at 
the start and end times of a flare. 

GOES observations of a B7 flare which occurred 1986 
January 15 at 10:09 UT, are shown in Figured] In the 
first column, the flare signal is assumed to dominate, 
i.e. background is set to zero, while in the second col- 
umn, the flare signal has been extracted by subtracting 
the pre-flare flux from the original lightcurve. The top 

8 http : //hesperia. gsf c .nasa.gOv/rhessidatacenter/complementary_data/i 
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Fig. 3. — GOES lightcurvcs and associated temperature and emission measure profiles for a B7 flare which occurred on f986 January 
15. The profiles in Figures [3ji-[3]i are not background subtracted. The profiles in Figures [3]5-[3fi have had the pre-flare flux in each channel 
subtracted, while Figures l3"t-f3t show the profiles obtained using the TEBBS method. The error bars represent the uncertainty quantified 
via the range of background subtractions found acceptable by TEBBS. 



row (Figures [2b and [3]s) shows the non-background sub- 
tracted lightcurvcs with the background levels overplot- 
ted as horizontal lines. The second row (Figures [3^ and 
[3jf) shows the lightcurves after background subtraction. 
(N.B. since the background in the left column is zero, Fig- 
ures \3jp, and[3f> are the same.) The third and fourth rows 
show the temperature and emission measure profiles, re- 
spectively, derived from the lightcurves shown in the sec- 
ond row. An acceptable temperature profile is shown in 
Figure [31c which peaks at 8 MK around 10:13 UT. How- 
ever the corresponding emission measure (Figure de- 
creases at the time of the flare. The reason for this un- 
characteristic behavior is that the background flux com- 
ponent is dominating the emission measure evolution of 
the flare, thereby making it impossible for properties to 
be derived accurately. Conversely, by subtracting the 
pre-flare flux, as shown in Figures [3]e-[3pi, significant arti- 
facts are introduced to both the temperature and emis- 
sion measure profiles. This is because this background 
subtraction is causing the flux ratio at the beginning of 
the flare to be comprised of two small numbers, which 
leads to large discontinuities when folded through the 
temperature and emission measure calculations. 

A more accurate approach would be to assume that 
the flare flux may also contain some contribution from 
the quiescent plasma from which it originates as shown 
in Figure [2] This assumption was the ba sis for the 
backg round subtraction method developed bv lBornmannl 
(1990). This technique applies three tests to a given com- 



bination of long and short channel background values: 
the increasing temperature test, the increasing emission 
measure test (together known as the increasing prop- 
erty tests), and the hot flare test; to determine whether 
a given choice of background levels produces physically 
meaningful results. The increasing property tests assume 
that both temperature and emission measure exhibit a 
characteristic overall increase during the rise phase. In 
these tests, background levels were selected and a pre- 
liminary subtraction was made. The relationship be- 
tween the long and short channel fluxes during the rise 
phase was approximated with a linear fit of the form, 
Fs = mF^ + c, where m is the slope and c is the in- 
tercept. From these fitted values, the temperature and 
emission measure for each point along the rise phase were 
were c alculated using the polynomials of iThomas et al.l 
(1985) and compared to their previous value. If over- 
all increases in these parameters were observed, then the 
background subtraction was said to have passed the in- 
creasing property tests. 

To pass the hot flare test of iBornmarml (|1990l ) , the 
background temperature (calculated by plugging the ra- 
tio of the background val ues, Rb = F^/F^i m t° t ne 
temperature polynomial of lThomas et al.lll985f) must be 
less than the background subtracted flare temperature at 
all times during the flare. This helps prevent unphysi- 
cal temperatures/emission measures being derived if the 
short channel approaches t he de tection threshold. 

The tests of IBornmarml £l990) were the first attempt 
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to isolate a GOES flare signal from the background con- 
tributions based on the validity of the results produced. 
However, they have their d rawbacks. The tests use the 
simple parameterizations of iThomas et al.l (|1985l ) to cal- 
culate temp e ratur e and emission measure. Since then, 
IWhite et al.l ()2005D devised tables from updated detector 
responses from GOES-l to GOES-12 and plasma source 
functions that take into account the marked differences 
in temperature and emission measure when derived us- 
ing co ronal rather than pho tospheric abundances. Previ- 
ously, IThomas et al.l (|1985[) had only provided one set of 
coefficients to their parameterization. In addition, Born- 
mann's tests do not take into account the GOES instru- 
mental temperature threshold. This threshold stands at 
4 MK and exists because such a temperature would cor- 
respond to a flux ratio of R = 1/100 which is beyond the 
sensitivity of the XRSs used in this study. This means 
that these tests may not always identify the background 
combinations which may lead to unphysical profiles. An- 
other shortcoming lies in the linear fit to the rise phase 
used in the incre asing proper t y test s. When demonstrat- 
ing the method, IBornmannl ( 19901) did not include the 
beginning of the rise phase in the linear fit because sig- 
nificant flux increases are often not observed there (e.g., 
FigureUJa) and can affect the fit's accuracy. However, this 
leaves the beginning of the rise phase untested, which is 
the period most likely to exhibit spikes or discontinuities 
due to an unsuitable background subtraction. If these 
spikes are big enough, they can easily be mistaken for the 
true peaks and produce unreliable results. In the next 
section the TEBBS method is described in detail and 
the ways in which it improves upon the above-mentioned 
short-comings of the Bornmann tests are discussed. 

3.2. Temperature and Emission measure-Based 
Background Subtraction (TEBBS) 

TEBBS has been developed to facilitate accurate cal- 
culation of the plasma properties of large numbers of 
solar flares observed by GOES. This is done by automat- 
ically isolating the flare signals from their background 
contributions. Bornmann's method has been updated 
and improved in a number of ways. Firstly, explicit 
tem perature and emiss ion measure values calculated us- 
ing IWhite et al.l (2005) are utilized in the background 
tests. This is favored over simply using the flux ra- 
tios because the characteristic temperature and emission 
measure evolution of the flare can be directly analyzed. 
Secondly, extra criteria have been added to the hot flare 
test so that the minimum background subtracted flare 
temperature must be greater than the instrumental tem- 
perature threshold of 4 MK. Similarly the maximum 
background subtracted temperature must be less than 
the upper limit of the IWhite et al.1 ()2005h tables, i.e., 
<100 MK (corresponding to a flux ratio <1). This up- 
per limit is much higher than any GOES temperatures 
found by previous studies. This helps to identify all pos- 
sible flare signals which produce unphysical profiles in- 
cluding those with discontinuities. Finally, another cri- 
terion has been added to the increasing property tests 
requiring that any temperature/emission measure value 
taken from the early rise phase which is not used in the 
linear fit, must be less than the peak taken from the 
rest of the rise phase. This helps to remove possible flare 
signals which show spikes at the beginning of the temper- 
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Fig. 4.— GOES XRS lightcurves from 1986 January 15 06:35- 
10:55 UT. The start and end times of the B7 flare shown in Fig- 
ures \3\ and [5] as defined by the G OES event list are marked by the 
dashed and dot-dashed vertical lines respectively. 



ature / emission measure profiles which may not be iden- 
tified by the original Bornmann tests. Both the TEBBS 
method and the original Bornmann tests assume that 
the background level in each channel is constant during 
the flare. This may not necessarily be the case, espe- 
cially when the flare occurs during the decay of an earlier 
event. This was deemed to be a rare enough occurrence 
(236 out of 34,361 events between 1991 and 2007, i.e., 
0.7% - see Section |2~2"T) that it would not introduce any 
significant errors. Moreover, as the peak flux and peak 
temperature occur near the beginning of an event, the 
slope of the background would have a negligible effect. 

The assumption underlying TEBBS is that the bound- 
ary between flare flux and backgrou nd lies som e where 
between zero and the pre-flare flux. IBornmannl ()1990f ) 
justified this by assuming a quiescent flux component 
from the flare plasma. However, there are a number of 
reasons why the background may not be well represented 
by the pre-flare flux. For example, if the recorded start 
time of the flare is later than the true start time, the 
flare flux will have already risen considerably, thereby 
causing the pre-flare flux to be much higher than the ac- 
tual background. Furthermore, if the flare occurs on the 
decay phase of another flare or at a time of high back- 
ground flux, the flare will not be seen in the XRS data 
until the flare flux dominates the flux from the rest of the 
solar disk. Thus the flux at the reported flare start time 
(i.e., pre-flare flux) is a convolution of background and 
early flare flux, and therefore should not be subtracted 
in its entirety. This was the case for the B7 flare on 1986 
January 15. GOES lightcurves from an extended period 
around the flare (06:45-10:55 UT) are shown in Figure|4j 
The start and end times as defined by the GOES event 
list of the B7 flare are shown as the vertical dashed and 
dot-dashed vertical lines respectively. It can be seen that 
this flare has occurred on the decay phase of an M-class 
flare which began around 06:50 UT. Because of this high 
pre-flare flux, the initial evolution of the B7 flare was 
not readily detectable in the XRS data. Therefore the 
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Fig. 5. — Short channel flux versus the long channel flux for the 
1986 January 15 flare (solid curve). The grey shaded area in the 
bottom left hand corner represents the possible combinations of 
background values from each channel for this event. The orange 
line represents a linear least-squares fit to the rise phase of the 
event. 

pre-flare flux was not an accurate approximation of the 
background and thus only a certain fraction should be 
subtracted. 

In order to apply the TEBBS method to this or any 
other flare, the first step is to define a sample space of 
possible background combinations. The range in each 
channel is between zero (equivalent to no background 
subtraction) and the minimum flux measured during the 
flare (equivalent to subtracting the pre-flare flux). The 
grey region in the bottom left of Figure [5] shows the sam- 
ple background space for the 1986 January 15 flare. This 
sample space is divided into twenty equally linearly sep- 
arated discrete values in each channel (F^ , ^J 3 ), thereby 
creating four hundred possible background combinations. 
Results were found to be independent of this binning and 
so twenty was chosen minimize computational time while 
ensuring that the background space was adequately sam- 
pled. Each background combination is then subtracted, 
thus creating four hundred sets of background subtracted 
lightcurves as possibilities for the flare signal. It is to 
these lightcurves that the hot flare test and increasing 
property tests are applied. 

The first test to be applied is the hot flare test. The 
minimum temperature, T m j n , of each lightcurve is cal- 
culated. Any background combinations corresponding 
to temperature profiles with a minimum temperature of 
Tmin < 4 MK are discarded. Then T m i n is compared 
with the background temperature, Tg, calculated using 
the background values, (F^,Fg). If T min < T B then 
that background combination is discarded. Furthermore, 
should the flux ratio at any point be greater than or equal 
to unity (i.e., T > 100 MK) the background combination 
is also discarded. The background combinations of the 
1986 January 15 flare which passed (solid region) and 
failed (hashed region) the hot flare test are shown in Fig- 
ure [B^. From this panel it can be seen that the number 
of possible background combinations has already been 
halved. 
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Fig. 6. — Sample background space for 1986 January 15 flare. 
The black shaded areas illustrate the range of values which pass 
a given background test, while the hashed regions denote back- 
ground values which fail, a) the hot flare test, b) the increasing 
temperature test; c) the increasing emission measure test, and d) 
points which passed all three, or failed one or more. 

Next, the increasing property tests are applied. As in 
the Bornmann tests, the relationship between the short 
and long channel fluxes during the rise phase is fitted 
with linear function of the form, Fg = mF^ + c, so as to 
reduce the influences of fluctuations in the data. Such a 
fit is justified by the fact that 90% of flares in this study 
have a Pearson correlation coefficient greater than 0.85 
for their rise phases and 95% of flares hav e a correla- 
tion c oefficient greater than 0.75. Following Bornmann 
(|1990l) . the first sixth of the rise phase duration was not 
included in the linear fit (orange line, Figure [5]). This is 
because significant increases are often not observed di- 
rectly after the GOES event list start time (see Figure [lj 
which can affect the accuracy of the fit to the rise phase. 
The choice not to include the first sixth was determined 
by experiment in order to exclude any non-increasing 
portion of the early rise phase. Using these linearly 
approximated values, the evolution of the temperature 
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Fig. 7. — Temperature and emission measure profiles for the 1986 January 15 flare for all possible background combinations. The left 
column shows profiles which passed all three tests, while the right column shows profiles which failed one or more tests. 



and emission measure during the rise phase is calculated 
from all of the background subtracted lightcurves. Each 
value is compared with its preceding one and the per- 
centage of times when the temperature/emission measure 
increases is calculated. Background combinations which 
result in profiles exhibiting a total increase less than a 
certain threshold are discarded. This threshold was cho- 
sen heuristically to be the maximum rise percentage from 
all four hundred background combinations minus seven 
e.g. if the maximum number of times that an increase 
in temperature was observed as a percentage of the to- 
tal rise phase is 77%, the threshold would be 70%. If 
this threshold leaves no background combinations which 
pass all three background tests it is iteratively reduced in 
steps of five percent until there is at least one background 
combination which passes all three tests. This method 
still leaves the very beginning of the rise phase untested. 
Therefore the next step is to calculate the temperature 
and emission measure profiles for the whole rise phase 
from the non-fitted background subtracted lightcurves. 
Any profiles that show a peak in the 'untested' section of 
the rise phase greater than the peak found in the 'tested' 
section are discarded. Figure [Hfr and [6]c show the back- 
ground combinations which passed (solid regions) and 
failed (hashed regions) the increasing temperature and 
increasing emission measure tests respectively. 

Having completed these tests, only background combi- 
nations which pass all three are deemed suitable. This 
leaves a small distribution of allowed background com- 
binations, shown as the solid region in Figure |SJi. The 
temperature and emission measure profiles correspond- 
ing to each of these background combinations are shown 
in Figure [7] (smoothed for illustrative purposes). The 
left column shows the profiles corresponding to back- 
ground combinations which passed all three tests, while 
the right column shows profiles corresponding to combi- 



nations which failed one or more tests. It can be seen that 
all the profiles in the left column are well behaved and 
are more conducive to calculating peak values and peak 
times. In contrast, many of the profiles in the right col- 
umn exhibit discontinuities and spikes, particularly near 
the beginning of the flare. It is impossible to calculate 
useful peak values and times from these profiles, either 
automatically, or even manually. Although some of the 
temperature profiles in the right column appear to be 
well-behaved, the corresponding emission measure pro- 
files contain artifacts, and vice versa. 

The TEBBS method was applied to the 1986 January 
15 B7 flare and the resulting time profiles are shown in 
the third column of Figure[3]([3j-[3^). The background lev- 
els were chosen from the combination closest to the center 
of the pass distribution in Figure |6H and are shown as the 
horizontal dashed and dot-dashed lines in Figure O. The 
background subtracted fluxes can be seen in Figure [3J/. 
The error bars mark the uncertainty in the background 
subtraction which was taken as the range of the pass 
distribution in Figure |6H. The TEBBS temperature and 
emission measure profiles are shown in Figures [3]fe and[3j, 
respectively. Both of these profiles show smooth rise and 
decay phases. Note that the temperature profile does 
not have a discontinuity as in Figure \5jg. Furthermore, 
the emission measure evolution is no longer dominated 
by the background contribution as in Figure EH nor does 
it exhibit spikes or discontinuities as in Figure GJfi. The 
error bars in these panels represent the range of accept- 
able temperature and emission measure profiles seen in 
the left column of Figure [7] Note that the uncertainties 
at the beginning of the flare are largest. This is expected 
as the flux ratio during the early rise phase is made of 
smaller numbers than the rest of the flare. Therefore, 
a slight inaccuracy in the background subtraction can 
cause a more significant change in temperature and emis- 
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sion measure. The beginning of the flare also shows an 
unexpectedly high temperature of >6 MK. This can be 
explained by the fact that the start time defined by the 
GOES event list was probably after the actual start time 
due to the high emission from the M-class flare which 
preceded it (see Figure QJ . This would imply that the 
flaring plasma was initially cooler than 6 MK but by the 
time the flare emission could be detected over that of the 
M-class flare, the plasma had already been heated sub- 
stantially. Figure [3] shows that TEBBS has performed a 
successful, automatic background subtraction, superior 
to either of those performed with the other two methods 
discussed in Section |3~T1 

Having successfully tested TEBBS on other flares cho- 
sen at random, the method was applied to all 51,196 
selected flares in the GOES event catalog from 1980 
January 1 to 2007 December 31. The specific back- 
ground combinations were chosen in the same way as 
for the 1986 January 15 flare. Of these, successful back- 
ground subtractions could not be performed for 1,870 
events (~3%) and so were discarded from the dataset. 
Of these 144 were 'double flares' or even 'triple flares' 
making the assumptions of TEBBS invalid. The remain- 
ing 1,726 can be characterised in the following ways: 
events dominated by 'bad' points; events where the short 
channel approaches the lower detection threshold leading 
to unphysical flux ratios and hence derived properties; 
unphysical light curves resembling square or triangular 
waves with up to an order of magnitude flux amplitude; 
events which we would be interested to analyze. The 
size distribution of events discarded by TEBBS was very 
similar to that of the original dataset. This shows that 
TEBBS did not preferentially discard flares of a partic- 
ular size and therefore did not bias the results of this 
study. The associated plasma properties (peak tempera- 
ture, peak emission measure, radiative loss rates, and to- 
tal radiative losses) were derived for the remaining 50,703 
events. Uncertainties on the plasma properties for each 
event were calculated from the corresponding range of 
allowed TEBBS background subtractions as was done in 
Figures |3]f-[3j. Values from 'bad' data points and neigh- 
boring data points were ignored due to their tendency 
to produce unreliable spikes when folded through the 
temperature and emission measure calculations. 'Bad' 
points are marked as such by the GOES software as un- 
trusted measurements because of the instrument states 
reported by telemetry (e.g., because of gain changes) or 
because they are outliers from surrounding data. The 
identification of these 'bad' points is justified by simul- 
taneous observations from other GOES spacecraft. The 
TEBBS database is publicly available on Solar Monitoi0 
and will also be available through the Heliophysics Inte- 
grated ObservatorjQ- The statistical relationships be- 
tween the above derived properties are discussed in the 
next section. 

4. RESULTS 

Peak temperature, peak emission measure, and total 
radiative losses each as a function of peak long channel 
flux are shown as a density of points in Figure [8] Each 
column shows distributions obtained using each of the 

9 http : //www . solarmonitor . org/TEBBS/ 

10 http://www.helio-vo.eu/ 



three background subtraction methods discussed in Sec- 
tion [3l Uncertainties on each data point are not shown 
for clarity. 

The relationship between peak flare temperature and 
peak long channel flux is shown in Figures [HMDs. While 
the non-background subtracted distribution in Figure Eh 
displays some trend of larger flares exhibiting higher tem- 
peratures, there is a flattening of the distribution below 
the CI level. This is due to the influence of the back- 
ground contributions which become highly significant at 
low fluxes. In addition, there is a stark vertical edge at 
Bl level due to the absence of A-class flares in the GOES 
event list. A horizontal edge at ~5 MK is also seen, 
which is due to the instrumental detection limit. There 
is more scatter in the pre-flare background subtracted 
distribution in Figured]!), with events of all classes show- 
ing temperatures in excess of 25 MK up to a temperature 
of ~80 MK (beyond the range of the plot axis) . By sub- 
tracting all of the pre-flare flux, the value of the flux ratio 
at the beginning of the flare can become erroneously large 
due to dividing one small number by another. This can 
lead to spuriously high temperature values when folded 
through the temperature calculations and can often be 
greater than the real peak temperature (see Section [3]). 
Many of the high temperature values (>25 MK) in Fig- 
ure particularly those corresponding to low peak long 
channel fluxes, have been taken from such spikes early 
in the flare. In contrast, the TEBBS method performs 
background subtractions which tend not to cause such 
temperature spikes. As a result the distribution in Fig- 
ure (5]c shows much less scatter. The flattening at low 
fluxes seen in Figure HJa has also been reduced by the use 
of the TEBBS method. 

In order to examine the relationship between peak tem- 
perature and peak long channel flux, a number of meth- 
ods were used. First, the Kendall tau coefficient was 
calculated which is a non-parametric correlation coef- 
ficient i.e. it does not assume a pre-defined model for 
the data. It is based on the rank of the data points 
rather than the values themselves, making it more suit- 
able than other correlation coefficients (e.g. the Pearson 
linear coefficient) to distributions with significant out- 
liers or scatter such as those in this study. The Kendal 
tau correlation coefficient for Figure was found to be 
0.42 which represents a statistically significant correla- 
tion. Next, the relationship was quantified using linear 
regression. Ordinary least squares (OLS) was not used 
however, because of three characteristics of the data: the 
presence of several outliers which produced non-normal 
behavior between an OLS regression fit and the observa- 
tions (i.e. the residuals were not normally distributed); 
data truncations due observational cutoffs below Bl level 
and 4 MK; the underlying power-law number distribution 
of the observations, i.e. the greater number of smaller 
events relative to larger ones. To address these char- 
acteristics of the observations, the methods of robust 
statistics were used. The basic assumption in OLS is 
that the residuals are normally distributed and the solu- 
tion to the problem is calculated by minimizing the sum 
of the squared residuals. However in this case, the sum 
is replaced by the median of the squared residuals. This 
results in an estimator that is resistant to the outliers 
by fin ding the narrowes t strip covering half the observa- 
tions (jRoussee uw 1984). To account for the population 
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Fig. 8. — Two-dimensional histograms of peak temperature, peak emission measure, and total radiative losses, each as a function of peak 
long channel flux, derived using various background subtraction techniques for all selected GOES events observed between 1980 and 2007. 
The data in the left column has not had the background subtracted. The data in the middle column had the pre-flare flux subtracted, while 
the distrib utions in the right column have been derived u sing the TEBB S metho d. Overplotted on p anels c and / ar e the r elationships 
derived by Garcia & Mcintosh (1992, purple dotted line), Feldman ct al. (1996b, orange dotted line), Battaglia ct al. (2005, red dashed 
line), and Hannah ct al. (2008, green dot-dashed line). The solid black lines show relationships described by Equations [5] E] and [8] The 
arrow heads mark events which are upper or lower limits due to saturation of the XRS and point in the direction that the true value would 
have been located. The crosses mark events which display lower flux limits and derived properties which are neither upper nor lower limits 
but only rough estimates. This is because their derived properties are functions of the fluxes in both channels which each saturated. See 
Appendix for a more detailed discussion of the effects of saturation on derived flare plasma properties and the saturation levels of the XRS 
on the various GOES satellites. 



distribution, the regression analysis was weighted using 
the flux values themselves, with smaller events weighted 
less than larger events. Th e truncation in the data i s 
handled using the method of Bhattacharya et al. (1983). 
The form of the fit resulting from this method is given 

by 

T = a + piog w F L MK (5) 

This form was chosen because it implies a linear relation- 
ship between temperature and the logio o f the peak long 
channe l flux such as those found b y both Feldma n et al.l 
(|1996bf) and lBattaglia et al.l ([20051) . The values of a and 
$ were found to be 34±3 and 3.9±0.5 respectively. Fi- 
nally, the goodness of this fit was examined by using a 
modified, robust R 2 statistic which quantifies the vari- 
ance in the data explained by the model. Whereas the 
usual R 2 value is based on the mean-squared-error, the 
modified robust R 2 statistic is based on the median (con- 
sistent with the robust fitting method used above). It 
also accounts for the degrees of freedom and in the fit- 
ting and the uncertainties on each data point. It was 
found that the modified robust R 2 value for the above fit 
was 0.62. This value is lowered by the structure in the 
distribution at least in part caused by the instrumental 
truncations below B-class and 4 MK. Nonetheless this 
value still implies that the Equation [5] is a suitable fit to 



the distribution. 

The relationship between peak emission measure and 
peak long channel flux is shown as a density of points 
in Figures (5]ff-[S]f. The non-background subtracted dis- 
tribution in Figure [5]cZ displays the same selection effect 
as that in Figures [SJa, i.e., a sharp cut-off at low GOES 
class (~B1 level). This cut-off is not seen as clearly in 
Figures [Us or [8]f) due to background subtraction. Large 
amounts of scatter were found below the Ml level in Fig- 
ures (5]i and [5]s which is not seen to the same degree 
in the TEBBS distribution in Figure The unusually 
high emission measures in Figures [8]i and [8]e have been 
recorded from erroneous features such as those in Fig- 
ures U3i and [3pi. There is a well defined linear lower edge 
in all thr ee of the dist ributi ons. A similar feature wa s 
found by [Garcia (1988]) and IGarcia" & Mclntosbl (1992). 
This edge is a natural consequence of the way emission 
measure is calculated, approximated by Equation [2j The 
second term in this equation asymptotically tends to zero 
which means it varies very little at high temperatures. 
As a result, emission measure becomes directly propor- 
tional to long channel flux causing the observed linear 
edge, which correspo nds to high temperatures. This fea- 
ture was also seen by IGarcia" fc Mcintosh] (|1992f ) but not 
explained. 

To examine the correlation between peak emission 
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Fig. 9. — Two-dimensional histograms of peak emission measure and total radiative losses as a function of peak temperature for each of 
the three background subtraction techniques. The arrow heads mark events which are upper or lower limits due to saturation of the XRS 
and point in the direction that the true value would have been located. The triangles mark events whose values are only rough estimates 
due to both channels saturating. See Appendix for a more detailed discussion of the effects of saturation on derived flare plasma properties 
and the saturation levels of the XRS on the various GOES satellites. 



measure and peak long channel flux, the Kendall tau 
coefficient of the TEBBS distribution in Figure |8]f was 
calculated and found to be 0.8, implying a significant 
correlation. In order to compare our results with those 
of previous studies, a linear relationship between the 
login of these two proper ti es, such as tho s e fou nd by 
Garcia fc Mcintosh! (11992ft . iBattaglia etaLl (|2005f) . and 
Hannah et al.l ( 2008ft . was applied to the data. The fit 
was performed in log-log space using the same linear re- 
gression method used for the relationship between the 
temperature peak and the long channel flux. To remain 
consistent with previous studies this fit was re-expressed 
as a power-law of the form: 



EM = 10 7 F? cm 



(6) 



The values of 7 and S were found to be 53±0.1 and 
0.86±0.02 respectively. This relationship is expressed in 
the inverse as 

F L = r)EM c W in" 2 (7) 

where rj and C were found to be 1 x 10- 61±1 and 1.15±0.02 
respectively. The modified robust R 2 value for the above 
model was found to be 0.73, implying a good fit. 

Total radiative losses as a function of peak long chan- 
nel flux is shown as a density of points in Figures \Sjg [5J. 
All three distributions clearly show an increasing trend 
with peak long channel flux. The similarity between the 
three distributions suggests that total radiative losses are 



not as sensitive to background subtraction as either peak 
temperature or peak emission measure. This is to be 
expected since peak values are taken from single points 
which can be very sensitive to erroneous spikes caused by 
inappropriate treatment of the background. However, to- 
tal radiative losses are integrated over the flare duration. 
Therefore, if a flare contains erroneous temperature or 
emission measure spikes, their contribution to the total 
radiative losses will not be as significant if the rest of the 
flare is 'well-behaved'. For small flares, however, the ef- 
fect of these erroneous values would be expected to have 
a greater influence. The 'turn-up' at A- and B-class lev- 
els in Figure Hfi is consistent with this. This distribution 
was found to have a high Kendall tau correlation coef- 
ficient of 0.73. It was then fit using the same method 
as used for the emission measure peak long channel flux 
relationship. The resulting fit was expressed in the form: 



Lrad = 10 £ F£ ergs 



(8) 



This form was chosen because it implies a linear rela- 
tionship between the logio of the total radiative losses 
and long channel peak flux with an intercept of e and a 
slope of 9. This is justified by a high Pearson correlation 
coefficient, calculated in log-log space as 0.8. The values 
for e and 9 were found to be 34± 0.4 and 0.9±0.07 re- 
spectively. The modified robust R 2 statistic was found 
to be 0.71, implying that Equation [8] well represents the 
distribution. 
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Distributions of peak emission measure and total radia- 
tive losses as a function of peak temperature are shown in 
Figure [9l Each column corresponds to distributions ob- 
tained from the same background subtraction methods 
as in Figure |5] 

Peak emission measure as a function of peak temper- 
ature is shown as a density of points in Figures [HMUc. 
A clear relationship between the two properties is not 
apparent in the non-background subtracted distribution 
in Figure EJa. A horizontal edge from 5-12 MK and just 
above 10 49 cm~ 3 is exhibited with the majority of flares 
located just below this edge. Any relationship between 
these two properties is even less clear in Figure 0>- Very 
large scatter extends beyond the range of this plot to 
~80 MK. The artifacts introduced into both the temper- 
ature and emission measure profiles by each of the respec- 
tive background subtraction methods (such as those in 
Figures [3jg andjU'i) have exacerbated the scatter. A more 
discernible trend with less scatter is revealed by the use 
of TEBBS in Figure [9]c. This distribution clearly shows 
that flares with hotter peak temperatures have greater 
peak emission measures. However, there seems to be an 
edge to this distribution at low temperatures which may 
also be explained by a limit of Equation [2] 

Total radiative losses as function of peak temperature 
is displayed as a density of points in Figures IHtHSF- Al- 
though the TEBBS distribution in Figures [9]f appears 
comparable to the non-background subtracted distribu- 
tion in Figure [9]i, it displays less scatter than seen in 
the pre-flare subtracted distribution in Figure |9]e which 
has data points extending beyond the range of the plot 
axis to ~80 MK. No clear trend between peak temper- 
ature and total radiative losses is discernible in any of 
Figures [9]i-[9]f, although Figures [9]rf and[9]f do show a ten- 
dency for higher temperature flares to have greater total 
radiative losses. This implies there is no strong relation- 
ship between these properties. This is despite the fact 
that total radiative losses are function of both temper- 
ature and emission measure. This absence of a trend 
may be due to the assumptions used in deriving these 
properties (e.g. constant density). 

5. DISCUSSION 

The TEBBS distributions in Figures [5] and [S] consis- 
tently show the least scatter and most discernible trends 
between properties derived from GOES data. The non- 
background subtracted and pre-flare subtracted distri- 
butions show a higher number of artifacts such as edges 
and anomalously high values. This shows that TEBBS 
is a superior method of automatically subtracting back- 
ground than either of the other two methods; first be- 
cause it successfully separates the flare signal from the 
background contributions, and second, produces fewer 
artifacts in doing so. However, there still may be biases 
in the distributions derived using TEBBS. Such biases 
may be due to the fact that TEBBS uses full-disk in- 
tegrated observations. A comparison between temper- 
ature and emission measure profiles produced in this 
study and those derived from spatially resolved instru- 
ments could further highlight how reliable the TEBBS 
results are and be used to quantify any systematic bi- 
ases. Spatially resolved observations could be taken 
from instruments which observe in similar wavelength 
bands to the XRS, such as the Soft X-ray Telescope 



(SXT) onboard Yohkoh, the Soft X-ray Imager (SXI) on- 
board GOES-12 and GOES-13, or the X-Ray Telescope 
(XRT) onboard Hinode. Such a study would be useful 
in further determining the strengths and weaknesses of 
the TEBBS method. Furthermore, it must be acknowl- 
edged that several necessary assumptions were used in 
calculating the plasma properties, as there are any time 
these properties are derived usin g GOES observations . 
In this study, coronal abundances (jFeldman et afl ll992'). 
a constant d ensity of 10 10 cm" 3 , the ionization equi- 
librium from iMazzotta et al.l (j!998f) , and an isothermal 
plasma were assumed. It has been shown that coro- 
nal iron abundances du ring flares can reach eight times 
the photospheric l e vel (jFeldman et al.lf2004f l. or higher, 
and IPhillipseTaTI (poTot T found coronal densities above 
10 13 cm" 3 using high-temperature density sensitive ra- 
tios. Using either of these assumptions in the calculation 
of the flaring plasma properties could affect the results. 
However, this was not done in this study so that these 
results would be more comparable with those of previous 
studies. 

The dist ribution in Figu r e [8]c w as co mpared with the 
studie s of iFeldman et al.l (|1996b| ) and feattagli a et al.1 
(2005). These studies found a linear correlation between 
peak temperature and the logio of peak long channel flux. 
In both papers, the relation was expressed in the form: 

F L = 3.5 x 1()P T +* W in" 2 (9) 

The values of /3 and n from these studies can be found 
in Table [1] Values from this study were calculated by 
rearranging Equation [5] into the form of Equation [§] 
and are also included in T able [T] for comp a rison. The 
IFeldman etall (|1996b[ ) and IBattaglia etall (|2005f ) rela- 
tions are also overplotted on Figure [8]c as the orange dot- 
ted and red short dashed lines respectively. The distri- 
bution of this study reveals predominantly lower temper- 
atures for a given long channel flux than both previous 
studies . Ther e is closer agreement with IFeldman et al.l 
(|1996bf ) than IBattaglia etall (|2005h for B-class events 
but beyond this, lowe r temperatures are obta ined. This 
can be explained by IFeldman et al.1 (|1996 b) using the 
BCS to calculate temperature. In that paper, it was 
stated that temperatures obtained with the BCS agreed 
with those from GOES below 12 MK but above this 
point were higher on average by a factor of 1.4. To in- 
vestigate this, the mean peak temperature of all flares 
in our sample of M-class or greater was computed and 
found to be 16 MK. This flux threshold was chosen 
for two reasons. Firstly, the difference between back- 
ground s ubtracted (this study) and non-background sub- 
tracted (jFeldman et al.l[l996b| ) results are negligible in 
this regime. Secondly, of these events, 95% had peak 
temperatures greater than 12 MK. This was compared to 
the mean Feldman temperature obtained for these events 
by plugging their long channel peak flux into the fit of 
IFeldman et al.l (|1996bf) . The Feldman mean temperature 
was found to be 20.9 MK which differs from that of this 
study by a factor of 1.3. This i s lower than the value 
quoted bv IFeldman et al.l (|1996bl ). This is because they 
measured temperature at the the time of the long channel 
peak which would be expected to occur after the temper- 
ature peak. Assuming that the temperature peaks before 
the long channel flux, the difference in ratios would im- 
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TABLE 1 

Values for Fl-T relationship (Equation [9} 



Study 



Fcldman ct al. (1996b) 
Battaglia et al. (20051 
TEBBS 



0.185 -9 
0.33±0.29 -12 
-9+. 1 



26+ - 04 
u - zo -o.oi 



TABLE 2 

Values for linear edge in Fl-EM distribution (Same form 
as Equation [6j 



Study 


7 


8 


Garcia & Mcintosh (1992) 


53.04 


0.83 


TEBBS (Eqn|B) 


53±0.1 


0.86±0.02 


TEBBS (lower edge) 


53.4 


0.96 



ply that a flare's temperature (M-class or greater) drops 
by 10% before the long channel peak. This implication 
is supported by this study, as part of which the tem- 
perature at the time of the long channel peak was also 
calculated. This shows that between the temperature 
and long channel peaks, a flare's temperature drops on 
average by 10%-11% for flares greater than or equal to 

M-class and 9%-10% for all flares. 

The slope of the relation of iBattaglia et al.l (|2005l ) ap- 
pears closer to that of the fit found by this study. How- 
ever, the relation consistently gives temperatures which 
are 3-4 MK higher. Th is disc repancy in the intercept 
is due to IBattaglia et~aH (|2005Q 's use of RHESSI to ob- 
tain the temperatures. The value of T in that study was 
calculated as either the temperature of an isothermal fit 
or the lower of two temperatures in a multi-thermal fit 
to a RHESSI spectrum. This was compared to GOES 
temperature, Tq, and a relation was derived given by 

T = 1.12T<3 + 3.12 MK (10) 

Substituting this into Equation |9] and rearranging into 
the form of Equation El values for a and (3 are found to 
be 28l; 1 5 8 and 2 - 7 ^i 7 3 3 respectively. Although these val- 
ues are similar to those found for this study (a = 33 ± 3; 
/3 = 3.9 ± 0.5), the large uncertainties mean that little 
statistical significance can be assigned to this similar- 
ity. This hi ghlights that a more comprehensive study 
than that of IBattaglia et al.l (|2005l) was needed to more 
precisely understand the statistical relationships between 
the thermal properties of solar flares. 

Next, the emission m easure distribution i n Figu re [8]f 
was compared work of iGarcia fc Mcinto sh (1992j). In 
that paper, the linear edge to this distribution was ad- 
dressed (also discussed in Section 2] of this paper) and a 
fit t o this lower b ound was quoted from a previous pa- 
per, IGarcial JT988). This was of the same form as Equa- 
tion [6] Garcia's values are shown in Table [2] and the fit 
corresponding to them is overplotted on Figure |8]f as the 
purple dotted line. However, this relation does not fit 
the lower bound of this distribution very well. In fact it 
seems to better fit the distribution itself being as it is so 
similar to the values found for Equation [6] in Section [4] 
(shown in the second row of Table [2]). A rough fit to this 
edge shows it is much better formulated by the parame- 
ters shown in the third row o f Table El The discrepancy 
may be because the sample of lGarcia fc Mcintosh! (|1992f ) 
was insufficient to reveal the actual limit of this relation- 
ship. However, it may be also be due to the fact that 
methods different from those of iWhite et al.l (|2005l ) were 
used to calculate tem perature and emission measure (e.g. 
iThomas et "all 119851 ). The credibility of this limit is im- 
portant as it suggests a well defined minimum amount of 
material emitting in the GOES passbands produced by 
a flare of a certain long channel peak flux. Although this 
limit is due to the nature of Equation [2] the result should 



TABLE 3 

Values for EM-F l relationship (Equation 



Study 




n 


c 


Battaglia ct al. 


(2005) 


3.6xl0" 50 


0.92±0.09 


Hannah ct al. i 


:2008) 


1.15X10" 52 


0.96 


TEBBS 




lxl0~ 61±1 


1.15±0.02 



be compared with statistical studies using other instru- 
ments to confirm whether it is a breakdown in the validity 
of the temperature and emission measure calculations of 
IWhite et al.l (pool or has any physical significance. 

This distribu ti on w as als o compared to t he w ork of 
IBattaglia et~"aH (|2005l ) and lHannah et all (120081) who 
found correlations between RHESSI emission measure 
and background subtracted GOES long channel peak 
flux. These relations were expressed in the same form as 
Equation [7] The values found by these studies are dis- 
played in Table [3] along with the values from this study 
for comparison. These previous fits are also overplot- 
ted on Figure |8]f as the dashed red and dot-dashed green 
lines respectively. These fits are steeper than our distri- 
bution. The relation of lHannah et al.l (|2008l ) however, 
agrees well at B- and C-class which was the range on 
which that study focused (A to low C-class). The differ- 
ence in slope can not be d ue to the different sen sitivities 
of GOES and RHESSI as lHannah et all ((2008T) showed 
that GOES emission measure is consistently a factor of 
two greater than that obtained from RHESSI. The differ- 
ence in slope may therefore be attributed to the extension 
of the distribution to M- and X-class. H owever, it may 
also h ave been affected by the fact that lHannah et al.l 
( 2008) calculated the emission measure at the time of the 
peak in the RHESSI 6-12 keV passband rather than the 
pe ak emission me a sure, as in this study. The relation 
of Battag lia et al.1 (120051 ) gives consistently lower emis- 
sion meas ures than both the T EBBS distribution and 
rel ation of lHannah et all (|2008T >. This can be explained 
bv IBattaglia et al.l (120051 ) measuring the emission mea- 
sure at the time of the hardest HXR peak which tends 
to occur early in the flare before the SXR and emission 
measure peaks. 

The distribution of peak emission measure as a 
function of peak temperature in Figure |H]c shows 
that hotter fla r es have larger peak emission measures. 
iFeldman etaLl (11996b! ) found a power-law relationship 
between these two properties expressed by 

EM = 1.7 x io 013T + 46 cm- 3 (11) 

However, the Pearson correlation coefficient between 
temperature and the log of emission measure in Figure[S]c 
was calculated to be 0.3. This implies that these proper- 
ties are in fact not well linearly correlated at all despite 
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an apparent tr end of hotter f l ares h aving higher emis- 
sion measures. Hannab Tet al.l (|2008l ) also examined the 
relationship between emission measure and temperature 
for A-low C-class flares and found no correlation. If the 
Figure [9]c distribution is examined more closely, there 
does not appear to be any r elationship between t he two 
properties within the range lHannah et all (|2008t ) stud- 
ied. Indeed, the Pearson correlation coeffi cient for C1.0- 
class and below is only 0.2. This supports lHannah et al.l 
(2008) findings. However further examination of this re- 
lationship is necessary to draw firmer conclusions. 

6. CONCLUSIONS & FUTURE WORK 

A method, TEBBS, has been presented for isolating 
the solar flare signal from GOES soft X-ray lightcurves 
by accounting for contributions from both solar and non- 
solar backgrounds. This allows the properties of the flar- 
ing plasma itself to be more accurately derived. It can be 
systematically applied to any number of flares, removing 
many of the inconsistencies that can be introduced when 
manually defining a background level. This makes it a 
particularly suitable method for conducting large-scale 
statistical studies of solar flares characteristics. TEBBS 
was found to produce fewer spurious artifacts in the de- 
rived temperature and emission measure profiles for both 
individual events (Figure |3]) and in large statistical sam- 
ples (Figures [8] and E|) , compared to when either all or 
none of the pre-flare flux was removed. This led to more 
reliable relationships being derived between flare plasma 
properties (temperature, emission measure etc.), which 
can in turn place constraints on the 'allowed' values of 
properties for a flare of a given GOES magnitude. 

TEBBS was successfully applied to 50,056 flares from 
B-class to X-class, making it the largest study of the 
thermal properties of solar flares to date. It was found 
that peak temperature scales logarithmically with peak 
long channel flux as described by Equation [5j Mean- 
while, peak emission measure and total radiative losses 
scaled with peak long channel flux as power-laws given 
by Equations [6] and [H Uncertainties were calculated for 
these derived relati ons unlike previous st udies. The ex- 
ception to this was Battagli a et al.l (|2005l ) who provided 
uncertainties for their slopes only. The uncertainties de- 
ri ved using TEBBS wer e nonetheless smaller than those 
of lBattaglia et al.l (|2005| ) and include uncertainties on the 
intercepts as well as slopes. Furthermore, while these re- 
sults are broadly in line with previous studies, it was 
found that flares of a given GOES class have lower tem- 
peratures and higher peak emission measures than pre- 
viously reported. 

Peak emission measure and total radiative losses were 
also examined as a function of peak temperature. It was 
found that flares with high peak temperatures al so have 
high peak emission measures (in ag reement with iGarcial 
119881 and iGarcia k. Mcintosh! 119921) . However, the de- 
rived correlation was relatively weak. Similarly, it was 
also found that flares of a given peak temperature could 
exhibit a large range of radiative losses with no clearly 
defined trend. This lack of a clearly defined relationship 
between two derived properties could be attributed to the 



assumptions that go in to calculating them. Although 
both a constant density and a fixed coronal abundance 
were assumed in this study, b oth have been shown to vary 
during individual events (e.g. lGraham et aT1l2011| ). A fol- 
lowup analysis of how changes in these variables might 
affect the derived properties, particularly in conjunction 
with hydrodynamic simulations, may lead to more rea- 
sonable correlations. 

This compilation of solar flare properties represents a 
valuable resource from which to conduct future large- 
scale s tatistical studies of fl are plasma properties. For ex- 
ample, IStoiseZeFjjl] (|2008| ) derived analytical predictions 
of temperature and emission measure in response to elec- 
tron beam and conduction driven heating and compared 
the results to RHESSI observations of 18 microflares. 
They found an order of magnitude discrepancy between 
conduction driven emission measures predicted by the 
Rosner- Tucker- Vaiana (RTV; iRosner et al.lll978T ) scaling 
laws and observation. This seemed to suggest that elec- 
tron beam processes dominated. However, they noted 
that RHESSFs high temperature sensitivity (>10 MK) 
mean that the observed temperatures may not have well 
represented the conduction value of the microflares, thus 
explaining the discrepancy. The fact that the GOES XRS 
has a sensitivity to lower temperatures than RHESSI 
makes the TEBBS database ideal for exploring this pos- 
sibility. Since the RTV scaling laws and electron beam 
heating models are widely used to understand and model 
solar flares, it is important to examine disagreements be- 
tween their predictions and observation. 

Another example of the use of RTV s c aling laws in 
understanding flares is lAschwanden et all (|2008|) . They 
used these laws to derive theoretical (EM oc T 4 3 ) 
and observed (EM oc T 4 7 ) scaling laws between peak 
temperature and emission measure for solar and stel- 
lar flares. However, as p art of their study, re sults 
from previous s t udies such as lFeldman et al.l (jl996b) and 
iFeldman et al.l ()1995|) were included which did not ac- 
count for background issues. TEBBS can therefore also 
be used to examine these scaling laws with greater sta- 
tistical certainty and therefore provide more clarity on 
the discrepancies between theory and observation. As 
the scaling laws derived by lAschwanden et al.l (|2008t ) ap- 
ply to solar and stellar flares, conclusions drawn from 
TEBBS can be extended to stellar flares as well. 

TEBBS can be used to examine a wide range of flare 
characteristics, such as thermodynamic evolution and, in 
light of the work of IStoiser et al.1 ()2008h . even flare loop 
topologies. As TEBBS is also the largest database of 
thermal flare plasma properties to date, it will provide a 
valuable resource for future solar flare research. 
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TABLE 4 
GOES Saturation Levels 



GOES Satellite 


Time Period 


Long 


Channel Flux 


Short Channel Flux 






(10 


-4 W m -2) 


(1CT 4 W m" 2 ) 


GOES-6 


06-Nov-1980 - 17-Dec-1982 




(...) 


1.8 


GOES-6 


24- Apr- 1984 




13 


1.2 


GOES-6 


20-May-1984 - 24-Jun-1988 




(...) 


1.2 


GOES-6 


06-Mar-1989 - 02-Nov-1992 




12 


1.2 


GOES-10 


02-Apr-2001 - 15-Apr-2001 




18 


4.7 


GOES-12 


28-Oct-2003 - 7-Scp-2005 




17 


4.9 



APPENDIX 
GOES SATURATION LEVELS 

During the period of this study, there were 32 X-class flares in the GOES event list which saturated either the short 
channel or both channels. No events included in this study saturated the long channel without also saturating the 
short channel. Saturation of the GOES channels has different and important effects when deriving each of the flare 
plasma properties. If the short channel saturates but the long channel does not, then the derived temperature during 
the period of saturation is a lower limit because T oc F$/ ' Fl, to a first order approximation. However, derived emission 
measure during the same period is an upper limit because EM oc T~ x , to a first order approximation. Likewise, since 
dL ra d/dt oc EM, radiative loss rates (and total radiative losses) are also upper limits. If both channels saturate 
however, then it cannot be determined (without extrapolation of the lightcurves) whether properties derived during 
the saturation period are upper or lower limits since they are all functions of the flux ratio. 

Within the GOES event list used in this study, only the XRSs onboard GOES-6, GOES-10, and GOES-12 were seen 
to saturate. Each channel in each XRS had different saturation levels which can be seen in Table 2] These values 
were taken from the GOES lightcurves of saturated events throughout the GOES event list. GOES-6 had a very long 
lifetime and as a result, the saturation levels of each channel were seen to degrade over time, also shown in Table |H 
(Dates are inclusive.) 
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